Universality for the breakup of invariant tori in Hamiltonian flows 
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In this article, we describe a new renormalization-group 
scheme for analyzing the breakup of invariant tori for Hamil- 
tonian systems with two degrees of freedom. The transfor- 
mation, which acts on Hamiltonians that are quadratic in the 
action variables, combines a rescaling of phase space and a 
partial elimination of irrelevant (non-resonant) frequencies. 
It is implemented numerically for the case applying to golden 
invariant tori. We find a nontrivial fixed point and compute 
the corresponding scaling and critical indices. If one compares 
flows to maps in the canonical way, our results are consistent 
with existing data on the breakup of golden invariant circles 
for area-preserving maps. 

PACS numbers: 05.45. -Kb, 64.60. Ak 



I. INTRODUCTION 

In Hamiltonian systems with two (or more) degrees 
of freedom, smooth invariant tori typically persist under 
small perturbations [|l]-Q . The most stable tori appear to 
be the ones for which the frequency ratio is a quadratic 
irrational, such as the golden mean. As the strength of 
the perturbation passes some critical value, these tori are 
observed to exhibit self-similar scaling behavior, and then 
they break up |^,^. To study such critical behavior, sev- 
eral renormalization schemes have been proposed [[/[-p^ . 
The idea is to use a transformation Tl that maps a Hamil- 
tonian H into a rescaled Hamiltonian TZ{H), in such a 
way that irrelevant degrees of freedom are contracted. 
The transformation TZ should have roughly the following 
properties: TZ has an attractive integrable fixed point Hq 
that has a smooth invariant torus of a given frequency. 
Every Hamiltonian in its domain of attraction 2? has a 
smooth invariant torus of the desired frequency. Another 
non-integrable fixed point lies on the boundary of 
V, also called critical surface, and this nontrivial fixed 
point attracts the Hamiltonians on the critical surface. 
Figure 1 shows schematically the expected nature of the 
renormalization flow in the space of Hamiltonians. The 
existence of a nontrivial fixed point has strong implica- 
tions concerning universal properties associated with the 
breakup of invariant tori. The analysis of the renormal- 
ization for area-preserving maps gives support to the 
validity of this general picture. 

The renormalization we define, following the scheme 
proposed in Ref. pd|], is similar in spirit to the block 
spin transformation in statistical mechanics, in the sense 
that it uses a process of "elimination and rescaling". 
The elimination of irrelevant frequencies is done by us- 



ing canonical transformations as in Kolmogorov-Arnold- 
Moser (KAM) theory. The frequencies we want to elim- 
inate are the non-slow modes (non-resonant part of the 
perturbation), i. e. the modes which only affect the mo- 
tion for a short time. The slow modes (resonant part 
of the perturbation), which produce the small denomi- 
nators in KAM theory, are shifted towards the non-slow 
modes by a rescaling. They can thus be eliminated by 
iteration. The rescaling is the same as in Refs. [pT|-p^. 
It includes a shift of the resonances, and a rescaling of 
the actions and of the energy. Here and in what fol- 
lows, the word "resonances" refers to the frequency vec- 
tors Vk = {Pk,q_k) defined by the continued fraction ap- 
proximants Pk/qk for the frequency ratio of the invariant 
torus. The corresponding closed orbits accumulate at the 
invariant torus |14|, and are often used in numerical in- 
vestigations. 

We consider the following class of Hamiltonians with 
two degrees of freedom, quadratic in the action variables 
A = {Ai,A2) and described by three scalar functions of 
the angles ip — {t^i, </'2) G (the two-dimensional torus 
parametrized by [0, 27r)^): 
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H{A,^)^-m{v){n-Ay 



(1.1) 



where is the frequency vector of the considered torus, 
and r2 = (1, a) is some other constant vector, not paral- 
lel to u>o- 

This family of Hamiltonians has been investigated in 
Refs. 15 1 . In particular, a KAM theorem was 

proven for this family in Ref. |l^ ] based on Thirring's ap- 
proach ]T^, in which the KAM transformations are con- 
structed such that the iteration stays within the space 
of Hamiltonians quadratic in the actions. In order to 
prove the existence of a torus with frequency vector 
for Hamiltonian systems described by Eq. (1.1), it is not 
necessary to eliminate m, but only g and /: One has to 
find a canonical transformation such that the equations 
of motion expressed in the new coordinates show trivially 
the existence of this torus. For Hamiltonians (1.1), if one 



takes g and / equal to zero, then the resulting equations 
of motion are 



dA 
~dt 



1 dm 

2 9^ 



{Q. ■ A) 



dtp 
'dt 



■m{(p){n ■ A)n + ujQ. 



(1.2) 



Then A = defines an invariant torus, and the motion 
on this torus is quasiperiodic with frequency vector ljq 
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(even if the resulting Hamiltonian is not integrable) . 
The elimination of / and g can be achieved with canon- 
ical transformations with generating functions that are 
linear in the acti on variables, and thus map the family 
of Hamiltonians (1.1) into itself. This is very convenient 
numerically, as one only works with three scalar func- 
tions m, g and /. The only approximation involved in 
our numerical implementation of the transformation is a 
truncation of the Fourier series of these functions, e. g. 
we approximate 



/(^) = 



by 



E 

veCL 



(1.3) 



(1.4) 



where Cl = {v & \\vi\ < L, ji^al < L}. 
In this paper, we focus on the frequency vector = 
(1/7,-1), associated with the golden mean 7 = (1 + 
\/5)/2. We choose the set of frequencies /~ describing 
the non-slow modes as the union of the two quadrants in 
the plane (I'l, 1^2) that contain the linear span of ujq: 



= {1/ e Z^li^ii^a < 0}. 



(1.5) 



The set /~ is depicted in Fig. 2. 

This "frequency cutoff" restricts to Fourier modes (the 
non-slow modes) that can be eliminated in one renor- 
malization step, without running into small denominator 
problems. As is common with cutoffs, there is not a sin- 
gle "natural" choice. Other possible choices of /~ will be 
mentioned later. 



II. KAM-RENORMALIZATION-GROUP 
TRANSFORMATION 

The renormalization scheme described in this section 
is for a torus of frequency vector ljq — (1/7,-1) where 
7 = {1 + \/5)/2. It is straightforward to adapt it to other 
reduced quadratic irrationals ||l7|| . 
Our transformation is composed of four steps: 

(1) A canonical change of coordinates, which acts on the 
Fourier mode e"' ''' for a resonance 1/ = i^'^, by a shift 
i/j, I— > i'k-2- The resonances for k < 2 will be eliminated 
in step 4, together with all the other frequency vectors 
in 

(2) A generalized canonical change of coordinates which 
corresponds to a rescaling of the action variables A. 

(3) A normalization corresponding to a rescaling of the 
energy (or time). 

The composition of the three steps described so far de- 
fines a map (m, g, /, a) 1-^ (m', g' , /', a'). 

(4) A KAM transformation (canonical change of vari- 
ables) that eliminates all non-slow modes from g' and 



/', i. e. the new functions /" and g" are in the range of 
I~, where I~ is a projection operator acting on a scalar 
function / of the angles as 



I" 



7(v) - E 



(2.1) 



Concerning step 4 we note that, for certain purposes it 
can be desirable to eliminate the non-slow modes from 
m' as well. But such a step generates terms of arbitrary 
order in the actions, which drastically complicates the 
analysis |l^ . 

We now give a more detailed description of the four steps 
that define the transformation. The first step is moti- 
vated by the observation that periodic orbits for the res- 
onant frequencies {ffc} accumulate geometrically at the 
(critical) WQ-torus. This suggests that the appropriate 
scaling (in the angle variables) is related to a shift in the 
sequence of resonances. In order to implement the shift 
i/fc 1-^ mentioned in step 1, we use the fact that the 
continued fraction approximants for the golden mean 7 
are 7^ = Fk+i/Fk, where Fk is the fc-th Fibonacci num- 
ber. These numbers can be defined recursively by the 
equation 



Fk+i 



= N 



Fk 
Fk-i 



with TV = 



1 1 
1 



, (2.2) 



starting with Fq = and Fi = 1. In other words, if we 
define Uq — (0, 1) then Vk — Nuk-i for A; > 1. Thus the 
desired shift Vk 1— > is induced by the linear trans- 

formation N~'^ , and the canonical transformation we are 
looking for is (A, ip) ^ {N~'^A, N'^^p). 
As the renormalization changes the scale, some of the 
slow modes become non-slow modes [e. g., (1, 1) is 
mapped into (0, 1) which is an element of /"]. We re- 
mark that there is no intrinsic sharp difference between 
the slow and non-slow modes. The boundary can be set 
at different places. We have chosen to include the coef- 
ficients (0, 1) and (1,0) among the non-slow modes (/~, 
to be eliminated by the KAM transformation), and (1,1) 
and (2, 1) among the slow modes. But conceptually and 
practically there would be no difficulty to include e. g. 
(1, 1) and (2, 1) (or any fixed finite number of resonances) 
among the non-slow modes. More generally, other choices 
in the splitting of {e^^'^} into slow and non-slow modes 
should lead to the same results provided e. g. that the 
ratio ji'l/lwo • ^\ is bounded on /~, and that con- 
tracts vectors v in the complement of /~ for some fixed 
A:>0|1|. 

The linear shift of the resonances multiplies Wo by 7 ^ 
(a;o is an eigenvector of iV); therefore we rescale the en- 
ergy by a factor 7^ in order to keep the frequency fixed 
at Wo- A consequence of the shift of the resonances is 
that n — (1,0;) is changed into Jl' = (1,q;'), where 
a' = {a + l)/{a + 2). 

Then we perform a rescaling of the action variables: we 
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change the Hamihonian H into 



HiA,<f)^XH[-,^ 



with A A(i/) such that the mean value of m' defined 
as 



(m') 



j2 {2n) 



(2.3) 



is equal to 1. Since the rescaling of energy and the shift 
iV^ transform the quadratic term of the Hamiltonian into 
7^(2 + afm{ip){n' ■ Af /2, this condition leads to A = 
7^(2 + af{m). 

In summary, the first three steps of our renormalization 
transformation rescale to, / and Q, — (1, a) into 



TO [^p) 



TO [N^'^^p) 



/'M=74(2 + a)2(TO)/(iV-V), 
, 1 + a 



2 + a' 



(2.4) 

(2.5) 
(2.6) 

(2.7) 



We remark that the map O i-^ O' given by Eq. (2.7) 
has an attractive fixed point fJ* = (1,7^^), which is an 
eigenvector of N"^ , with eigenvalue 7^ > 1. 
The fourth step is carried out via an iterative proce- 
dure, similar to a Newton algorithm. We start with 
Hq = H' . To simplify the description of the iteration 
step Hn 1-^ Hn+i, consider first the case where g„ and 
/„ depend on a (small) parameter e„ , in such a way that 
l~gn and I~/„ are of order ©(£„). The idea is to elimi- 
nate the non-slow modes of g„ and /„ to first order in e„, 
at the expense of adding terms that are of order 0{en) 
in the slow modes and of order 0(£^) in the non-slow 
modes. 

We will work with Lie 

transformations Un ■ (yn-i,A„_i) 1-^ ((^„,A„) gener- 
ated by a function Sn linear in the action variables, of 
the form 

Sn{A, if) = Yn{cp)n' ■ A + Zn{ip) + fln^' • , (2.8) 



characterized by two scalar functions y„, Z„, and a con- 
stant a„. The expression of the Hamiltonian in the new 
variables is obtained by the following equation : 

Hn+l = Hn oUn — c'^^" Hn 

= Hn + {Sn, Hn} + ^{Sn, {Sn, Hn}} + ■■■. (2.9) 

A consequence of the linearity of Sn in A is that the 
Hamiltonian Hn+i is again quadratic in the actions, and 
of the form 

Hn+i{A,ip) = ]^mn+i{p){Vl' ■ Af 



We notice that the vector Jl' remains unchanged from 
one step to the next. The functions Yn, Zn, and the con- 
stant ttn are chosen in such a way that gn+i and I^/n+i 
vanish to order 0(e„); see Appendix ^ for details. Con- 
sequently, I~gn+i and 1~ fn+i are of order ©(eq"), where 
Eq denotes the order of I^g' and I^/'. If Eq is small, 
an iteration of this procedure should define a canonical 
transformation Uh — U10U20 ■ ■ ■ oUn" ■ ■ ■ such that the 
Hamiltonian expressed in these new variables has only 
slow modes in g and /. In our numerical implementation 
this is indeed the case, even without small parameters. 
In summary, our renormalization-group (RG) transfor- 
mation acts as follows: First, some of the slow modes 
(resonant part of the perturbation) are turned into non- 
slow modes by a frequency shift and a rescaling. Then, 
a KAM-type iteration eliminates these non-slow (non- 
resonant) modes, while producing some new slow modes. 



III. DETERMINATION OF THE CRITICAL 
COUPLING; FIXED POINT OF THE KAM-RG 
TRANSFORMATION 

We start with the same initial Hamiltonian as in Refs. 

HiA, p)^^in-Af+uJo-A + efiif) , (3.1) 

where n = (1,0), loo = (1/7,-1), 7 = (1 + \/5)/2, and 
a perturbation 



f{<p) = COs(lVi • if) + C0S(Z^2 • f). 



(3.2) 



where v>i — (1,0) and v>2 = (1,1)- We represent all 
the functions by their Fourier series truncated by re- 
taining only the coefficients in the square Cl which con- 
tains {2L + 1)^ Fourier coefficients. For fixed L, we take 
successively larger couplings s and determine whether 
the KAM-RG iteration converges to a Hamiltonian with 
/ = 0, g = (trivial fixed point), or whether it diverges 
{f,g — > 00). By a bisection procedure, we determine 
the critical coupling ec{L). We repeat the calculation 
with larger numbers of Fourier coefficients, to obtain a 
more accurate approximation. Table I lists some values 
of £c{L)- For instance, for L = 10, the result is given 
with 7 digits; as a comparison, the method developed in 
Refs. ||l^,|l^ yields 4 significant digits (in order to obtain 
7 significant digits with this method one needs to cal- 
culate £c{L) up to L = 34). Moreover, we observe the 
disappearance of the oscillations of Fig. 1 of Ref. 
that is to say, the present method converges much faster 
than the one of Refs. [E2 13|. 



By iterating the transformation starting from a point 
on the critical surface, we observe that the process con- 
verges to a nontrivial fixed point H^ (or more generally 
to a nontrivial fixed set related to this nontrivial fixed 



[wo + g„+i(¥>)r2'] • A + fn+i{(p). (2.10) point by symmetries |13 2|]), which we characterize by 
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the Fourier coefficients of the three functions ^ g*, f* 
and r2* = (1,7^1). 

Figures 3 and 4 show the weight of the Fourier coefficients 
of m, . We observe that these coefficients decrease only 
slowly in the direction of u}q, and in particular, along 
the rescaled resonances N^^Vi, k > 0, which indicates 
that TO* is not analytic. This is due to the fact that our 
transformation does not eliminate the non-slow modes of 
TO (for the reasons mentioned earlier). 
Figures 5 and 6 show the weight of the Fourier coeffi- 
cients for the functions g^, and /», respectively. These 
coefficients seem to decrease exponentially, which indi- 
cates that and /* are real analytic. Notice also that, 
by construction, g^ and have only slow modes. 
At the fixed point iJ,, we compute the relevant eigen- 
values (critical exponents) for the linearized KAM-RG 
transformation. Since the breakup of invariant tori is 
observed in one-parameter families of Hamiltonians, in- 
dicating that the critical surface is of codimension one, 
one expects to find a simple eigenvalue S > 1, and no 
other spectrum outside the open unit disk. This is pre- 
cisely what we find numerically. Table I lists the values of 
(5 as a function of L. We obtain S G [2.650175, 2.650234], 
which is in agreement with the value obtained by MacKay 
for area-preserving maps ((5 = 2.650221). 
If the renormalization- group picture is correct, then all 
the relevant informations about critical tori with fre- 
quency vector u)o = (1/7,-1) is contained in the non- 
trivial fixed point H^. In particular, determines 
the observed critical scaling of phase space |^,^j8|. For 
the scaling factor A» — X{H^), we obtain numerically 
A, e [18.827910, 18.828203]; the value obtained for area- 
preserving maps is A* = 18.828171 (given in Refs. [^,^,^). 



IV. CONCLUSION 

We have shown that in the renormalization-group ap- 
proach to critical invariant tori, a partial elimination of 
non-resonant frequencies leads to substantial qualitative 
and quantitative improvements. Compared with previ- 
ous schemes that have been studied numerically, our new 
KAM-RG transformation yields more accurate results 
and a better defined domain of attraction for the non- 
trivial fixed point. Conceptually it is close to the type of 
transformations used for the study of critical phenomena 
in statistical mechanics. 

We have implemented our transformation for the case of 
tori with frequency vector u}q = (1/7,-1), where 7 is the 
golden mean. The extension to other reduced quadratic 
irrationals should be straightforward and yields similar 
results. It also seems possible to extend our transfor- 
mation to systems with three |2^ and more degrees of 
freedom [|ll| , where very little is known about the be- 
havior of invariant tori under perturbations that are not 
necessarily small 23 2^. 



ACKNOWLEDGMENTS 

We thank G. Bcnfatto, G. Gallavotti, R. S. MacKay, 
and P. Wittwer for helpful discussions. H. Koch acknowl- 
edges support by the National Science Foundation under 
Grant No. DMS-9705095. Support from the EC contract 
No. ERBCHRXCT94-0460 for the project "Stability and 
universality in classical mechanics" and from the Conseil 
Regional de Bourgogne is also acknowledged. 



APPENDIX A: ONE STEP OF THE KAM 
ITERATION 

In this appendix, we use the notation d f — df /dip. 
We consider a Hamiltonian of the form (1.1), with g — 
(1 - l-)G + et-G and / = (1 - l-)F + el-p. The 
parameter e is introduced for bookkeeping purposes only; 
it will be set to one at the end. Our goal is to find a 
canonical transformation U such that the functions /' 
and g' for the Hamiltonian H' = H o U have only slow 
modes, up to terms of order Oie^). We perform a Lie 
transformation defined in Sec. |ll[ The expression of t he 
Hamiltonian in the new variables is given by Eq. (^.9[): 
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l\2 



(Al) 



where 



g' = g + iVo-dY + m{n ■ dZ + ail^) + 0{e^), (A2) 
f = f + uJo- dZ + aujo-f^ + 0{e^), (A3) 

with Y, Z and a of order 0(e) to b e de termined. The 
generating function S given by Eq. ( ^.8| ) is determined 
by the projection of Eqs. (A2)-(A3) on the space of non- 
slow modes. We recall that the condition is that I^/' 
and I~g' are of order 0(2^). This leads to the following 
equations: 



I f + ljq ■ dZ = const, 

rg + ujQ-dY + r {mn ■ dZ) + l^man^ = 0. 
Equation ( [A5| ) determines a: 

(g) + {mn ■ dZ) 



a = 



(A4) 
(A5) 



(A6) 



Moreover the functions Y and Z have only non-slow 
modes, and are given by the following series: 



Wo • V 



-fue 



(A7) 



Y{^)= y ^—{g^ + {mn-dZ), 



(A8) 
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where /~* = I~ \ {0}. The transformed Hamiltonian 
(Al) is constructed by defining = H and i/'*-* for 
I = 1, 2, ... by the recursive relation 

ij('+i)(A,<p) = {5(A,¥p),i7W(A,cp)} 



which leads to 



+5 



H' 



+ (A9) 



E 

i=0 



(AlO) 



This can be expressed in terms of the image of the three 
scalar functions {m,g,f) given by the following equa- 
tions: 




(m(0),g(°),/(")) 



i'm,g,f), 



m 



(1) 



9 



(1) 



= 2TOr2 • ar - rr2 • dm, 

gft dY -Yfl- dg 



+mfl ■ dZ + maif + • dY, 
f^^^i = -YVl ■ df + gfl- dZ + gan'^ 
m(*+i) = 2m^'^n dY -Yn- dm^'\ 



LjQ ■ dZ, 



9 



9 



(*)rj dY -Yfl dg 



+mWrj • dZ + m^''>an^, 
for z > 1. 



(All) 

(A12) 
(A13) 

(A14) 
(A15) 
(A16) 

(A17) 
(A18) 



TABLE I. Parameters associated with the breakup of noble 
tori. 



L 


Critical value ec{L) 


Area multiplier 


Unstable eigenvalue 




0.027590 


18.828171 


2.650221 


3 


0.027625 


18.842654 


2.649660 


5 


0.027579 


18.823481 


2.652722 


8 


0.027588 


18.825034 


2.649653 


12 


0.027590 


18.829142 


2.650484 


13 


0.027590 


18.827464 


2.650082 


20 


0.027590 


18.828203 


2.650234 


21 


0.027590 


18.827772 


2.650151 


33 


0.027590 


18.828177 


2.650226 


34 


0.027590 


18.827910 


2.650175 



FIG . 1. Renormalization flow in the space of Hamiltonians 
(1.1) associated with the universal class of ujq = (I/7, —1). 
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FIG. 2. Non-slow modes (in the grey part) and slow modes 
(in the white part) associated with = (1/7' ^1)- 
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FIG. 3. Weight of the Fourier cocfRcicnts of m*: 
White: < 10-1°, grey levels: [10-i°, lO"'^], 
[10-^ 10-5], [10-^ 10-3], [10-^ 5.10-2], black: > 5.10-2. 
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FIG. 4. Weight of the Fourier coefficients of m*: 
M{ui) = max„2 |rn*y|. 
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FIG. 5. Weight of the Fourier coefficients of g*: 
White: < 10-1°, grey levels: [10-^°, lO-''], 
[10-^, 10-5], [10-^ 10-3], [10-3,5.10-2], black: > 5.10-2. 
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FIG. 6. Weight of the Fourier coefficients of /*: 
White: < 10-1°, grey levels: [10-^°, IQ-^], 
[10-^, 10-5], [10-5, 10-3], [10-3,5.10-2], black: > 5.10-2. 
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